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Trans-Planckian physics and signature change events in Bose gas hydrodynamics 
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We present an example of emergent spacetime as the hydrodynamic limit of a more fundamental 
microscopic theory. The low-energy, long-wavelength limit in our model is dominated by collective 
variables that generate an effective Lorentzian metric. This system naturally exhibits a microscopic 
mechanism allowing us to perform controlled signature change between Lorentzian and Rieman- 
nian geometries. We calculate the number of particles produced from a finite-duration Euclidean- 
signature event, where we take the position that to a good approximation the dynamics is dominated 
by the evolution of the linearized perturbations, as suggested by Calzetta and Hu [Phys. Rev. A 
68 (2003) 043625]. We adapt the ideas presented by Dray et al. [Gen. Rel. Grav. 23 (1991) 967], 
such that the field and its canonical momentum are continuous at the signature-change event. 
We investigate the interplay between the underlying microscopic structure and the emergent grav- 
itational field, focussing on its impact on particle production in the ultraviolet regime. In general, 
this can be thought of as the combination of trans-Planckian physics and signature-change physics. 
Further we investigate the possibility of using the proposed signature change event as an amplifier 
for analogue "cosmological particle production" in condensed matter experiments. 

PACS numbers: 04.62.+v; 98.80.-k; 47.37.+q; 67.57.De 



I. INTRODUCTION AND MOTIVATION 

Emergent spacetimes [J H, H| allow us to approach the 
subject of curved spacetime quantum field theory (CST- 
QFT) through the back door. This has been demon- 
strated in detail by using ultra-cold (non-relativistic), 
highly dilute Bose gases [3, IE El ■ Under appropriate 
conditions the fundamental microscopic theory can be 
replaced by a classical mean- field, the Bose-Einstein con- 
densate (BEC) 0. Collective excitations, both clas- 
sical and quantum, experience an effective spacetime 
whose entries are purely macroscopic mean-field vari- 
ables. The kinematic equations for linearized perturba- 
tions — neglecting back-reaction and finite temperature 
effects — are equivalent to covariant minimally coupled 
scalar fields, with a d'Alembertian 
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defined by an effective metric g a b 0, IE El • 
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Quantum field theory in curved spacetime is a good ap- 
proximation for semiclassical gravity — at the level where 
back-reactions of the quantum fields on the gravitational 
field arc negligible [H . As a consequence quantum effects 
in curved spacetimes (CSTs) — e.g., Hawking radiation 
[U EE EH and cosmological particle production — do not 
require emergent Einstein gravity per se; the existence 
of an emergent spacetime, an effective gravitational field 
g a b, is sufficient. Of course both systems, semiclassical 
gravity and any analogue model, must involve some dy- 
namics — and so they will eventually diverge from each 
other. The only possible loophole would be if Einstein 
gravity were itself to be the "hydrodynamics" of some 
more fundamental theory (of microscopic objects, for ex- 
ample strings, molecules, or atoms). For our purposes a 
perfect match is not required, and we refer the interested 
reader to 0, El Q EE EH for further details. 

It has been shown that the repulsive or attractive na- 
ture of atomic interaction in a Bose gas is directly re- 
lated to the signature of the low-energy emergent met- 
ric: Lorentzian (— , + + +) for repulsive interactions, Eu- 
clidean (+, + + +) for attractive interactions [3]. In 2001 
a BEC experiment [l?], EH was carried out that can be 
viewed as the first analogue model exp eriment. By tun- 
ing through a Fcshbach resonance [191 ] the atomic inter- 
actions were driven into a weakly attractive regime, and 
triggered a controlled condensate collapse. T wo y ears 
later, the theoretical work of Calzetta and Hu [2fj|. l2l| . 
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connected the the so-called Bosc-nova phenomenon with 
the amplification, mixing of positive and negative modes, 
and squeezing of vacuum fluctuations due to a signature 
change event. Perhaps surprisingly, the calculations car- 
ried out in 0, did not include background conden- 
sate dynamics and yet their theoretical predictions reflect 
the experimental data relatively well. For short time- 
intervals of attractive atom- atom interaction (i.e., a brief 
excursion into Euclidean signature) , the Bose-nova event 
is dominated by the evolution of the quantum perturba- 
tions, and to a good approximation independent of the 
background condensate dynamics. 

Specifically, given that we do not as yet have any pre- 
cise detailed model for emergent Einstein gravity, as op- 
posed to emergent curved spacetime, it seems necessary 
to focus on quantum effects that are merely of kinematic 
rather than dynamic nature. However, there are ways 
to study the influence of possible quantum gravity can- 
didates with CST-QFT. This branch of physics is called 
quantum gravity phenomenology (QGP) [221 ] - Emergent 
spacetimes can be used to analyze some portions of QGP, 
where Lorentz invariance violations (LIV) are present at 
ultraviolet scales. The LIV scale is supposed to be con- 
nected with the Planck length, where new physics is ex- 
pected. This idea is naturally implemented in any emer- 
gent spacetime model. For example, the BEC-based ana- 
logue models only recover Lorentz invariance (LI) for the 
low-energy, long-wavelength phonon modes. For trans- 
phononic modes microscopic corrections (e.g., quantum 
pressure effects) have to be absorbed into the macro- 
scopic picture. The borderline between the two modes 
can be viewed as the analogue LIV scale. Given that 
trans-phononic modes start to see first signs of the funda- 
mental microscopic theory, it makes sense to speak of LIV 
at the analogue Planck scale. The key reason why it is 
interesting to study modifications in the dispersion rela- 
tion is that many different effective field theories (EFTs) 
already predict deviations at the kinematical level. A de- 
tailed treatment of the analogue trans-Planckian model 
can be found in [H, H H, M, S3 • 

In ongoing work [28[ a classical phase-space method 
has been used to numerically simulate cosmological parti- 
cle production in BECs. There it is shown that for a con- 
sistent treatment microscopic corrections play an impor- 
tant role in the emergent spacetime picture. This leads 
to emergent "rainbow metrics" , with a Planck-suppressed 
momentum dependence for the modes, and consequently 
leads to a modification in the quasi-particle spectrum. 

Our central idea in the current article is to merge all of 
the above, to address the trans-Planckian problem for a 
signature change event in a Bose gas. We are particularly 
interested in the ultraviolet physics of the phonon modes 
and hence have chosen a specific BEC set-up where the 
external trapping potential does not interfere with the 
dynamics (e.g., hard- walled box). Step- by-step, in sec- 
tion (p"Tj) we show how spacetime emerges from a Bose 
gas; calculate the quasi-particle production from sudden 
sound speed variations in section (|III[) ; extend this calcu- 



lation to sudden variations for finite regions with differ- 
ent signature (Lorentzian <-» Euclidean) in section pVft : 
introduce ultraviolet physics and re-calculate the quasi- 
particle production in section ((Vj); suggest in section (|VI[) 
to employing a finite Euclidean region as a particle ampli- 
fier for cosmological particle production in a BEC; and 
last but certainly not least, we address the theoretical 
and experimental impact of our results in section (IVIII) . 



II. EMERGENT GEOMETRY FROM A BOSE 
GAS 

The intent in this section is to give readers unfamiliar 
with the topic some understanding of analogue space- 
times. Following the example of [l[ we introduce an ultra- 
cold (i.e., non-relativistic) , highly dilute and weakly in- 
teracting gas of Bosons, using the formalism of canonical 
quantization, and derive the equation of motion for small 
quantum fluctuations around some classical background; 
better known as a BEC, see section pi Al) . The Bosc- 
Einstein condensate is a state of matter where the Bosons 
macroscopically occupy the lowest quantum state. In 
section pi Bp we focus on the hydrodynamic case and 
recover a covariant minimally coupled massless scalar 
field for small quantum fluctuations in the BEC. Subse- 
quently, in section (fVj) . we will revisit this derivation and 
include trans-phononic modes (e.g., ultraviolet physics) 
into the emergent spacetime picture. 



A. Ultra-cold, weakly interacting Bose gas 

Suppose we have a system of N Bosons. In quan- 
tum field theory the field variables are quantum opera- 
tors that act on quantum states (Hilbert space of states; 
e.g., Fock space). Field operators either create, i/^^x), 
or destroy, ip(t,x), an individual Boson at a particular 
point in space and time, and satisfy the commutators: 



#,x),^,x') = ft{t,x),ft(t,x!) 



0; 



(2) 



#,x),Vit(t,x') =J(x-x'). (3) 



For a gas of trapped, ultra-cold, highly dilute and weakly 
interacting Bosons the Hamiltonian is given by 

H = Jdx (^,t|!_v 2 ^ + ftVrtd + j ftftw) ■ (4) 

This is a sum of the kinetic energy of the Boson field, and 
the two potential energy contributions; the external trap 
V cx t, and the particle interactions. The extreme dilution 
of the gas (e.g., 10 13 — 10 15 atoms/cm 3 ) suppresses more- 
than-two-particle interactions, and in the weakly inter- 
acting regime the actual inter-atom potential has been 
approximated by a pseudo-contact potential, 



U 



(5) 
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Here m is the single-Boson mass, and a the s-wave scat- 
tering length. The sign of the scattering length deter- 
mines the qualitative behavior of the interactions, 



a > 
a < 



repulsive ; 
attractive . 



(6) 



Negative and positive values of a are experimentally ac- 
cessible by tuning external magnetic fields, that inter- 
act with the inter-atomic potential; this process is called 
Feshbach resonance [l^] . We would like to emphasize the 
importance of Eq. (|6|) for the remaining sections. 

We now have all the necessary information about our 
system, encoded in Eqs. (J2J)— (|4j) , to calculate its dynam- 
ics. We use the Hciscnberg equation of motion to get the 
time-evolution for the field operator; 



.. dip 



(7) 



To apply this discussion to the emergent spacetime pro- 
gramme, we use the macroscopic occupation of the lowest 
quantum state below a critical temperature (e.g., for al- 
kali gases, below 10~ 5 K). If the cooling process prohibits 
the gas to solidify, a new state of matter will occur, the 
Bose-Einstein condensate. The condensate is a complex- 
valued macroscopic mean-field 

(ftt, x)) = ip(t, x) = V"oM exp (i<f> (t, x)) , (8) 

where the individual microscopic particles give way to 
collective variables. We will show that the condensate 
density, no, and phase, <po, define the analogue spacetime 
for small quantum fluctuations in the BEC. The essential 
step in deriving the kinematics for the perturbation is to 
separate the perturbation from the condensate: 



$ ~ V(*,x)+#(t,x); 
ft ~ ft(t,x) + 6ft(t,x). 



(9) 
(10) 



This transformation is canonical if the creation 6ft and 
destruction Sip of perturbations is consistent with the 
commutator Eqs. (J5J)- ©. Therefore, the commutators 
for the linearized quantum fluctuations are 



Sip(t,x),5ftt,x') = 6ft{t,x),5ft{t,x') =0; (11) 
Sftt,x),6ft(t,x')] =<5(x-x'). (12) 



The quantum perturbations are small perturbations 
which, as per Eq. ((5|), are represented by two collective 
parameters; the density n, and the phase <p. Clearly, any 
quantum perturbation should be related to variations in 
these two parameters; n ~ no + h and <p — 0o + </>■ A 
straightforward expansion of ijj and ft around uq and <po 
leads to 



Sip 
5ft 



1 2 n 



up 



2 no 



(13) 
(14) 



In this way Eqs. ([9]) and (JTUJ) are compatible with Eq. (8J. 
Thus the density, h, and phase, <f>, fluctuations operators 
arc Hcrmitian operators: 



(15) 



n 




' ' 8ft 




i 




( sj> 


~2 




{ 1> 



(16) 



It is easy to see that the new operators are a set of canon- 
ical variables: 



n(t,x),n(t,x')J = [0(i,x),0(i,x')J =0; (17) 
h{t, x), 4>(t, x')j = iS(x - x') . (18) 



The latter can be further modified, and we will subse- 
quently revisit Eq. (|T8|) . While the split into background 
plus perturbation is up to this point exact, we now lin- 
earize by assuming the perturbation to be small, allowing 
us to neglect quadratic and higher-order products of the 
perturbation field Sip. (There are also more sophisticated 
calculational techniques available based on the Hartree- 
Fock-Bogoliubov-Popov approximation, but they are an 
unnecessary complication in the present situation.) 

Merging equations l|9"j)-([TU]) with equations (|T5] ) -(fT3 1) . 
applying them to Eq. l[7|). we obtain two equations, 



dh „ 



n ° h Y7l , - 

V0 + n ■ v 

m 



0; 



^+V-V0+-71 = O. 

at h 



(19) 



(20) 



Here we introduce the background velocity of the con- 
densate, 



v = — V(/) - 
m 



(21) 



The quantity U can be thought of as an effective atomic 
interaction, as seen by the collective excitations. For 
long- wavelength, low-energetic modes this simplifies to 
Eq. (0, the usual pseudo-contact potential; 



phononic modes . 



(22) 



We will come back to this point in section (|V| where we 
shall write down the expression for U in general, and in 
the eikonal limit, in a manner appropriate to describe 
trans-phononic (or, ultraviolet) modes. In the language 
of condensed matter physics, we include "quantum pres- 
sure" effects. Without quantum pressure we restrict our 
analysis to the phononic regime, and assume that all col- 
lective excitations propagate with the same speed, the 
sound speed cq; 



n U 



(23) 
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Before we continue with our program, we would like 
to revisit the commutator derived in Eq. (fT8"]) . In the 
phononic regime, we are able to write Eq. (|20|) as follows, 



-uU +v -v*=-ui£- (24) 



Thus the commutator (|18p can be written in terms of <p 
and its fluid-following derivative (or material derivative) 
D(j>/Dt: 



where / is easily found to be, 



0(t,x),ft^(t,x')J =i<5(x-x') 
Here we have defined, 



TT -1EI 

* ~ U Dt 



(25) 



(26) 



which can be viewed as the conjugate momentum of 4>: 
see Eq. ([23]) below. 

The relation ([2"4"]) enables us to read density perturba- 
tions as fluid-following derivatives of the phase perturba- 
tions, and hence supplies us with the necessary tool to 
eliminate all occurrences of n from equation (|19[) : 



rab 



_ A 
f x 

u u v 

— V 



ngh _ h_ 2 

m h U 

~JfV x V z 



U u y 



~~V Vz 
-jjv x v z 

nph At, 2 



V 



v- 



The two equations, ([50]) . and 



are equivalent if 



(31) 



f ab = V|det( ffo6 )| g ab 



(32) 



Here g ab is a contravariant tensor, and since g ah g c b = 5 a c , 
it is only a question of matrix inversion to find its covari- 
ant equivalent, (l2"9"]h For considerably more details and a 
thorough derivation we suggest the following literature: 

a a am. 

At this stage we would like to comment on the physical 
implications of the results presented so far. The motiva- 
tion to write the differential Eq. (|2"T|) , governing the exci- 
tation spectrum in the form presented in (|28p . is to find 
an analogy for the curved spacetime Klein-Gordon equa- 
tion which describes minimally coupled spin-zero Bosons 
in curved spacetime. It is appropriate to define an emer- 
gent Lagrange density, 



-d. 





-ft 










-W0 


- V 





-V 



noh 



v<t> 



V V0 V 



(27) 



This equation, governing the kinematics for the phase 
perturbations <f>, is the connection between condensed 
matter physics and emergent quantum field theory in 
curved spacetimes. 



B. Analogue spacetime 

A compact and insightful way to express the evolution 
of phase perturbations (|27[) is 

1 



■s/\Aet(g ah )\ 
where we introduce 



d a [V\det(g ab )\g ab d b )^ = , (28) 



9ab 



CO 

u/h 



- [Co - v- 

—v x 



-v x 

1 






~Vy 



1 





-v z 



1 



(29) 

a covariant metric rank two tensor, whose entries are 
purely collective variables. The conformal factor depends 
on the spatial dimensionality, d, of the condensate cloud. 
To derive Eq. (|28|). we first write Eq. ([27]) as, 



d a (f ab d b $) 



0. 



(30) 



c = --f ab d a <i>d b < 



(33) 



such that Eq. ([28]) can be obtained as the Euler-Lagrange 
equations justified by the principle of least action. The 
momentum conjugate to cf> is specified by 



n; 



dC 

d(d t 4>) 



-f 



tb 



(34) 



and hence is in agreement with equation ([2"6"]) . We see 
that quantum phase and density perturbations in a Bose- 
Einstein condensate are a canonical set of field and con- 
jugate field operators on the emergent spacetime. 

We now temporarily set aside this analogy and apply 
our model to particle production in non-smooth emergent 
geometries. 



III. SUDDEN CHANGES IN SPACETIME 
GEOMETRY 

In this section we calculate the mixing of positive and 
negative frequencies due to "sudden" step-wise varia- 
tions in the sound speed, Eq. (|23[) . The initial and final 
emergent geometries are now flat Minkowski spacetimes, 
which are discontinuously connected at the step. Physi- 
cally the step is generated by a very rapid change in the 
magnetic field, which very rapidly drives one through a 
Fcshbach resonance, which in turn very rapidly changes 
the scattering length a, and so finally induces a rapid 
change in the speed of sound. 
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A. Quantum fields as harmonic oscillators 

An idealized Bose-Einstein condensate trapped in a 
finite quantization box of volume L d is comparable to flat 
Minkowski spacctime. Here the macroscopic parameters 
are zero background velocity, v = 0, and constant sound 
speed (|23p . That is, one considers a uniform number 
density, no{t, x) = const, and a fixed scattering length, 
a(t) = const. Therefore the emergent metric given in (j2"9")) 
is a diagonal tensor whose entries are time- and space 
independent; 



rab 



U 









n^h Q 
m 

Q n h 

m 

o o 



(35) 



We employ the canonical variables on our effective rel- 
ativistic spacctime to write the Klein-Gordon Eq. 

as 



dt [ JJ dt</> I — V </> = . 

U I m 



(36) 



It is possible to decouple the Klein-Gordon field into in- 
dependent Harmonic oscillators. To show this, we make 
use of 



</>(t, x) = 



1 



1 



L d/2 ^2 



E 



d k v* k (t)< 



ot«fc(*)< 



— ikx 



(37) 

Note that for a hard-walled box the modes fulfill non- 
periodic boundary conditions. In Minkowski spacctime 
there exists a natural set of mode functions, 



Vk(t) 



1 



(38) 



associated with the Poincare group, a symmetry group of 
the Minkowski line-element. Here d t is a time-translation 
Killing vector which can also be thought of as a differ- 
ential operator with eigenvalues (—iLu k ), where uo k > 
are said to be positive frequency modes. Hence the vac- 
uum is invariant under the action of the Poincare group 
and all observers agree on existence, or non-existence, 
of particles in flat spacctime. The physics is observer- 
independent, as expected. For a more detailed treatment 
see reference Q, and the appendix below. 

The decoupled equations for the mode-operators are 

a k {t) = -iu> k a k (t) ; and a' k (t) = +iuj k a\{t) , (39) 

where d k (t) = a k v k {t) and a\.{t) = a\v k (t). For now, we 
are working in the low-energy, long-wavelength regime 
and obtain a "relativistic" dispersion relation for the 
modes: 



Wk = ck, 



(40) 



Later on, in section {V}, we will include trans-phononic 
modes into our picture, and see how the microscopic 



structure induces LIV breaking terms for high-energy, 
short-wavelength perturbations. 

The canonical creation and destruction operators obey 
the usual commutator 



ak(t),a'_ k ,(t) =5, 



Jkk' 



(41) 



and acting on the particle basis for the Hilbert space of 
states, the Fock states | ), they arc a powerful tool to 
calculate the number of particles n k in the mode labeled 
by k. The normalized basis vectors can be obtained from 
the vacuum, or zero-particle state, |0). Thus is the state 
that is destroyed by operators fifc|0) = 0, for all modes 
labeled by k. In general 



and 



il\n k ) = {n + l) 1 ' 2 \{n + l) k ) 



ak\n k ) = n 1/2 \(n - l) fe ) 



(42) 



(43) 



so that we can define the number operator A/*,, 

N k \n k ) = a\a k \n k ) = n k \n k ) ■ (44) 

We now use the tools we have presented above to calcu- 
late the quasi-particle production in a BEC-based emer- 
gent spacetime due to sudden changes in the sound speed; 
such that we patch two flat spacetimes in a step- wise fash- 
ion, by suitably changing the microscopic parameters. 

B. Particle production 

From the emergent spacetime point of view step-wise 
changes in the microscopic parameters induce sudden 
variations in the collective variables, and enforce discon- 
tinuous interchanges between different Minkowski space- 
times. In terms of our emergent spacetime from a Bose 
gas, this can be achieved through an external magnetic 
field, that adjusts the atomic interactions, and conse- 
quently the scattering length a and the sound speeds c^. 

Following the ideas of Dray et al. [13, HH , we assume 
that the fields are continuous, 



^ =0. 

E 



(45) 



on the space like hypcrsurface J2 a t fixed time t = i*. 
Here t = is the time at which a transition occurs from 
& to (f. The Klein-Gordon equation in each region can 
be written in terms of exterior derivatives, 



{-l) n d a (T b d b ^j = *d{*d^) =: *dF i = ; (46) 

hence in each region dF- = 0, where F- is an exact (n— 1)- 
form. These forms are connected discontinuously; 



F = QF i +(l- Q)F± 



(47) 



6 



where is the usual Hcavyside function. Thus, we get 
dF = QdF i +(l-e)dFi-+S(t)dtA[F] (48) 



S(t) dt A [F] , 



(49) 



where [F] = F- — FK This supplies us with a connection 
condition for the canonical momentum, 



6{t) dt A [F] = ; 



IL 



Hi 1 



0. 



(50) 



(51) 



Thus the field and its canonical momentum must be con- 
tinuously connected across the spatial hypersurface of 
sudden variation. In the following we apply the con- 
nection conditions we have just found, and calculate 
the number of quasi-particlcs produced by sudden sound 
speed variations in our BEC. 

In each region i the mode expansion, 



i 



(52) 



involves a distinct set of creation ar k and destruction a\ 

i i * 

operators. The mode functions u~ k and u~ k obey the 
curved-spacetime Klein-Gordon equation, 



di 



dtu~ k (t) 



n h 7 2 



fcX(t)=0: 



(53) 



they form a complete basis for the two dimensional so- 
lution space. At the transition time 1M- we connect the 
mode functions and their first derivatives, 

Ui(t±) = M--Ui(&), (54) 
^d t Um = ^^#). (55) 

Here we have combined both mode functions into one 
mode vector U k , 



(56) 



In the cases we are interested in the transition matrix 
M 41 = M(0) is time-independent, and for a pair of com- 
plex conjugate mode functions in each region i and j, the 
matrix is of the form 



Q 



[ k Pk 



ft 



(57) 



The Wronskian W[u- k ,ir k ] = (d t ir k )ir k - u\ (d t ir k ) of 
the mode functions is time-independent, see Eq. (|53|) . 
which implies a constraint on the Bogliubov coefficients 



a k and (f k ; 



a 



-y.i 2 



(58) 



To start with, consider a no-particle state in the initial 
Minkowski spacetime (region i), such that ff^jO 1 ) = 0- 
The correlation between the mode operators can easily 
be constructed from (/^-(V 1 ) = (j^-(t^), considering Eq. (p>4"|) 
and Eq. (|55|) we get 



,t 



12 *3 

jj* ^j} oil* -i 

a k a k - Pk a - 



(59) 
(60) 



Consequently, the mean number of i-particles in the j- 
vacuum, 



(oipvi|oi> = |/#| 2 ^(o) 



(61) 



depends only on \£f k \. 

Explicitly, the elements for the transition matrix in a 
sudden change from region i to region j at a time t^- are 
given by: 




(62) 
(63) 



where to simplify formulae it is convenient to introduce 
the ratio of change in the dispersion relations X*, 



(64) 



The last expression in (|64|) is only valid in the hydro- 
dynamic approximation, where the number of particles 
produced from a single step is /c-independent. The mean 
number of particles in each mode fc is given by 



m 2 




(65) 



The advantage of this representation lies in the simple 
way that it can be extended for m sudden variations in 
a row 



M 12 . M M... Jlf (»-!)* 



(66) 



Mi 



and each single transition matrix M-^- +k ^ is of the form 
([57)) evaluated at t l<J - +1 \ The resulting matrix M-- car- 
ries the final Bogoliubov coefficients aqf and (3 k ~ for the 
whole chain of events. We would like to point out that 
the choice for the mode functions in the intermediate 
regimes does not have any influence on the final outcome, 
as can be seen in Eq. ([66]l . However for the validity of 
our calculation it is necessary to choose a pair of complex 
conjugate mode functions. The rest of the paper further 
investigates a particular scenario, that is two sudden vari- 
ations in a row, see Fig. [TJ In this situation 



a k 



a h a. 



■ft ft* 



12 nil r>12 2 

a k Pk +Pk a k 



(67) 

(68) 
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interactions 





a > 



a < 



1 1 * 



2 2 * 



3 3 * 



time 



FIG. 1: Two sudden variations in the strength of the atomic 
interactions is analogous to two space-like hyper-surfaces con- 
necting three different Minkowski spacetimes emerging form 
a Bose gas. 



To obtain the Bogoliubov coefficients for a two-step pro- 
cess, we use the single-step results (|fl2"j) and (|qU|) to ob- 
tain: 



07 



exp 



i cos(O fc T) - X + * sin(fi fc T) I , (69) 



and 



X 1 - 



cxp 
- 1 



cos(fi fc T) + sin(n k T) I ;(70) 



where for the finite duration intermediate region 2 we 
have defined the time interval T = t M — t—, and the dis- 
persion relation = u;7. The mean number of particles 
produced during this process is then given by 



13,2 (^-l) 2 



AX± 



(71) 



sin(ft fc T) 5 



It is easy to sec that the mean number of particles 
produced in such a process oscillates between two single 
step solutions of the form 



(x-i) 

AX 



(72) 



The upper bound replaces two up-down/ down-up steps 
with a single step (here given by X — X—/X-), and the 



lower bound replaces two up-down/ down- up steps with 
a single step (here given by X = X— X— ). We have illus- 
trated the quasi-particle production in the hydrodynamic 



limit obtained by two-sudden steps in a row in Fig. 4(a) 
In the following we show how to extend the calcula- 
tions for finite intervals during which the microscopic 
atoms experience negative (attractive) interactions. For 
small negative values of the scattering length a the con- 
densate description continues to exist for short periods 
of time, hence motivates the study of Euclidean emer- 
gent geometries from a Bose gas. A sign change in the 
atomic interactions corresponds to a signature change in 
general relativity. As a first step we investigate quasi- 
particle production from such an event in the hydrody- 
namic limit. 



IV. SIGNATURE CHANGE EVENTS 

In this section we investigate the behavior of emergent 
spacetimes arising from a Bose gas with variations in the 
principal nature of the interactions; for a finite amount 
of time we switch to attractive atomic interactions. For 
short time-scales and small absolute values of the attrac- 
tive s-wave potential, it is possible to hold on to the con- 
cept of an emergent spacetime. With the nature of the 
microscopic interactions switching from repulsive to at- 
tractive, the geometric hydrodynamics also changes from 
Lorentzian to Ricmannian. This supplies us with a toy 
model for signature change events, which we are going to 
investigate next. 

We first present the standard general relativity point 
of view, before adapting our previous calculations for the 
particle production due to sudden variations on finite du- 
ration Euclidean regions. 



A. Classical aspects 

Overall, we investigate manifolds that allow both Rie- 
mannian and Lorentzian regions. The latter is a gen- 
eralization of Minkowski spacetimes (special relativity), 
while Ricmannian geometries are Euclidean signature 
spaces including curvature (e.g., the surface of an or- 
ange). Distances in Riemannian spacetimes are positive 
semi-definite, while Lorentzian distances can be imagi- 
nary (time- like), zero (light-like), or positive (space- like) . 

In Minkowski spacetimes the signature can easily be 
read off as the sign of positive and negative eigenvalues 
of 1^^. From Eq. (|35[) we get Lorentzian signature 
( — , + + +). Arbitrary curved spacetimes are locally flat, 
and the signature can be read off from the pattern of 
eigenvalues of the metric tensor g ao at each point on the 
manifold. In Lorentzian signature the "time" coordinate 
can be chosen to have a different sign from the "spatial" 
coordinates. This is in contrast to the Riemannian, 
or Euclidean signature (+, + + +), where a distinction 
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between space and time as such does not exist [321 ] . 

We are mainly interested in the interface between these 
two spacetime geometries, and we investigate the physics 
around a space-like hypersurface ^ that separates the 
two spaces. There are two ways to be driven through a 
signature change, continuously or discontinuously. 

We are in favour of a non-smooth signature change, 
to avoid degeneracies of the effective gravitational field 
at the surfaces of separation. For continuous signature 
changes the metric volume element, and hence the ex- 
istence of an orthonormal frame vanish at J^; while for 
discontinuous signature changes the metric volume el- 
ement (hence, the orthonormal frame) is well behaved 

Thus we are allowed to transfer the connection condi- 
tions, as derived in our previous calculations for sudden 
changes in purely Lorentzian geometries (see Eq. (|54p . 
Eq. (|55j) , and the discussion in the appendix) , to now con- 
nect regions of Minkowskian and Euclidean spacctimcs 

M. 



B. 



Particle production from a finite-duration 
Euclidian region 



atomic 
interactions 



m = n 1 . 



a > 



a < 





1 






1 1 * 




u k' u k 







M- 



T 



3 3 * 




time 



FIG. 2: Two sudden variations in the strength of the atomic 
interactions are analogues to two space-like hyper-surfaces 
connecting three different spacetimes emerging from a Bose 
gas. While the initial and final geometries are Lorentzian, we 
allow the atomic interactions to be attractive in the interme- 
diate region, corresponding to an Euclidean geometry. 



In the following we repeat our calculations from section 
piip , but this time we choose an Euclidean geometry for 
region 2; see figure ([2]). As already pointed out above, 
the connection conditions remain the same. We further 
choose the mode functions in the Lorentzian regimes 1, 



and 3 to be as defined in Eq. (|38|) , where the sound speeds 
are well behaved; (cq) 2 ~ U > 0. For the intermediate 
Euclidean region we pick a special set of mode functions 
for (cq) 2 ~ U < 0, such that, 



1 



I 



(cosh(|fifc|i) +ismh(\Q. k \t)) ; (73) 
: (cosh(|fi fc |t) -isinh(|fi fc |t)) . (74) 



The two mode functions remain a complex conjugate 
pair. Here we make use of the purely imaginary dis- 
persion relation f2fc = i \ flk\- In general, with the trans- 
formation matrix S, 



S 



Vi((l-i) + 
2 + 



(75) 



we can map between the two different sets of mode func- 
tions. 



(76) 



Notice that S" 1 = S' t and det(S) = 1. 

It is easy to see that M— Uj; transforms to M— S Vr, 
while M— Ufc transforms to S^ 1 M^\j\. Altogether 



M— = M m S 5 _1 M M = M— M 1 



(77) 



is independent of the choice for the mode functions in the 
intermediate regime. 

Thus we are allowed to use our previous results, replac- 
ing nk -> i\flk\, and X— -> i \X^\ and X— -> -i \X% 
and keeping everything else. For the mean number of 
particles we now obtain 



\Pt\ 2 
1 

+ 4 



(X** - l) 2 



4X± 



\X 2 



\X' i 



(78) 

sinh(|^|T) 2 . 



Our results are compatible with Dray et al. [30 . l3l| , to 
the extent that the calculations and physical models over- 
lap, but we extend the calculation for arbitrary values of 
the sound speeds c~ . Furthermore we show that one can 
keep the calculations for purely sudden variation, as long 
as one picks a pair of complex conjugate mode functions 
in the intermediate regime. The standard general rel- 
ativity calculation corresponds to X— = i, X— = —i, 
X^ = l,so that 



M\ttk\Tf 



(79) 



The basic reason for this tremendous simplification is 
that in pure general relativity (with, by definition, a sin- 
gle unique spacetime metric) one always has the freedom 
to choose coordinates such that c = 1 in the Lorentzian 
region, and c = i in the Euclidean region. This is a free- 
dom we do not have in our BEC-based analogue space- 
time — the way this shows up in our calculations is that 
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a rapid change in the scattering length a has two ef- 
fects in the condensed Bose gas: First the speed of sound 
is changed, modifying the "signal cones"; and secondly 
the dimension-dependent conformal factor shifts by a fi- 
nite amount. This second effect is absent in the tradi- 
tional general relativity calculation of Dray et al. — for 
those authors it is sufficient to posit a specific and simple 
change in the metric tensor g a b and calculate the result- 
ing particle production. In our present situation, we first 
derive a specific (dimension-independent) change in the 
tensor density f ab induced by changing the scattering 
length, and then derive the corresponding (dimension- 
dependent) change in the metric tensor g a (,. The two 
situations are very closely related, but they are not quite 
identical. 

We have plotted the mean number of particles pro- 
duced in the Euclidean region in figure 4(c) The graph 



shows the quasi-particle spectrum as a function of k. 
The number of particles produced depends on fl k T. The 
longer the duration of Euclidian period, the more parti- 
cles will be produced during this process. 

However, there is a fundamental problem with the 
quasi-particle production in our effective spacctime. 
Given that it is possible to connect actual condensate 
excitations with the calculated quasi-particle spectrum, 
we expect the total number of particles produced to be 
finite. But in both cases, for sudden variations with and 
without signature changes, our results imply an infinite 
number for the total quasi-particle production; 



A. Rainbow geometries 

Calculations up to this point were based on the as- 
sumption that spatial variations in the overall conden- 
sate are small. More specifically, variations in the kinetic 
energy of the condensate are considered to be negligible, 
compared to the internal potential energy of the Bosons, 



h 2 V 2 y/n +n 



"2m y/no + 



(81) 



The left hand term in relation (|81[) is the quantum pres- 
sure term and is approximated to zero in the hydrody- 
namic limit . 

We now keep the quantum pressure term. A straight- 
forward computation shows that it is possible to absorb 
all first-order corrections into: 



U = U 



(Vn ) 2 - (V 2 n )n Vn 



V + V' 



(82) 

That is, the effect of quantum pressure can be absorbed 
into an effective atom- atom interaction, see Eq. (|22|). For 
a uniform condensate this further simplifies to 



U = U 



Amnr 



(83) 



In order to obtain an emergent gravitational field, we 
apply the eikonal approximation, 



N = 2 d - 1 n I a 



d-l 



(80) 



In the next section wc show how this problem resolves 
itself once microscopic corrections to the emergent space- 
time picture are taken into account. 



V ik 



u k = u + 



AmriQ 



(84) 



Thus also the speed of sound for ultraviolet modes has 
to be modified 



4 



U k n 



e qp k 



(85) 



Here we introduce, 



V. ULTRAVIOLET CORRECTIONS 

Up to now, we have restricted our calculations to the 
hydrodynamic limit, which is appropriate to describe 
the infrared behavior of the system. Low-energy exci- 
tations in the BEC are longitudinal phonon modes ap- 
proximately propagating with the same speed; see the 
dispersion relation given in Eq. (f40|) . It is well known 
that so-called trans-phononic modes show a non-linear 
relation between excitation energy and wavelength. In 
references [H, [H [H [H, H3 it has been pointed out 
that this kind of behavior might be viewed as ultraviolet 
corrections at the "analogue Planck scale" , the borderline 
between phononic and trans-phononic modes. 

Our next task is to use the eikonal approximation to 
include ultraviolet modes into the emergent spacetime 
picture. 



h 

2m 



(86) 



Wc define the borderline between the phononic and trans- 
phononic modes, 



1 e z 
x qp 

2 eg 



(87) 



to agree with the healing length £ of the condensate. The 
healing length is the distance over which localized per- 
turbations in the condensate tend to smooth out [35| . 
Phononic excitations have wavelengths that are much 
larger than the healing length, k S> £. These modes 
are relativistic modes in the sense of an emergent rel- 
ativistic dispersion relation (see Eq. (|4T)|) ). They prop- 
agate through an emergent gravitational field, given in 
(|2"9"|) . Higher energy excitations around the healing scale, 
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k ~ £, start to see deviations from the mean-field de- 
scription. Such modes do not only experience collective 
condensate variables, they start to see the bigger picture 
behind the mean-field. Consequently, they exhibit a non- 
relativistic dispersion relation, 



H) k = c k k 



elk 2 



£q p fc 4 . 



These microscopic corrections also influence the emer- 
gent spacetime, including an effective energy-dependent 
metric, 



9ab 



Ck 

u k /h 



(4 - v 2 

-v x 

~Vy 

-V z 



-V x 
1 






~Vy 



1 





-v z 



1 



which results in so-called rainbow geometries. 

Before we continue with our program, that is to re- 
calculate the mean number of particles produced includ- 
ing quantum pressure effects, we would like to empha- 
size that all the corrections are naturally small. For 
a detailed description of rainbow geometries, and the 
suppres sion mechanism of LIV terms in Bose gases see 

m IMS HzL S3- 



B. LIV and Particle production 

As mentioned above, the spacetime description only 
holds in the cikonal approximation. Nevertheless, a care- 
ful analysis for the particle production can be made, 
where the atomic interaction is a function of differen- 
tial operators, see Eq. (f52"| . For an isotropic condensate, 
the two different approaches will lead to the same set of 
equations, 



0, 



^ dt u\(t)y^k^ k ( t ) 



0. 



(90) 



The mode functions vr k are still formally represented by 
Eq. except that we replace to 1 with Q 1 , which now 

includes the e qp quantum pressure term. 

In the next two sections we connect those modes over 
three discontinuously patched spacetime geometries, and 
calculate the amplification and mixing of positive and 
negative modes living on such a manifold. 



1. LIV and L-L-L sudden 

We start with purely Lorentzian geometries, but allow 
two sudden variations in the speed of sound. The cal- 
culations are completely in analogy with that presented 
in section (|IVB[) . Here, the sound speeds are given by 
Eq. ([55]) , and from the definition for the ratio X* in the 
hydrodynamic approximation, we now define: 



as the sound-speed ratio for trans-phononic modes. 
Formally, the results for the Bogoliubov coefficients 
- a^r and /3j~, given in Eqs. ([69]) and (|70|) — and 
consequently the mean number of particles produced 
during this two-step process as given in Eq. (|7Tj) . all 
remain the same. The only effort required is to replace 
X^- with X^-, and urf, with ur k . 

W e have plotted the quasi-particle spectra in figure 
|4(b)| It can be seen that the quasi-particle spectrum for 
the two-step process is still oscillating between the two 
single-step processes represented by the black and green 
curve. Compared to the figure on the left, where we 
plotte d the same process in the hydrodynamic limit, see 
figure 4(a) we notice that the ultraviolet particle pro- 



(g9) duction rapidly approaches zero. 



2. LIV and L-E-L sudden 

Finally, we are left to analyze sudden variations from 
Lorentzian, to Euclidean, and back to Lorentzian space- 
times for rainbow geometries. It is interesting to notice, 
that in the presence of quantum pressure corrections the 
meaning of signature is also energy-dependent. We dis- 
tinguish between the following cases: 



1^1 < c o/ e qp 
Slg - 2 -^ \k\>4/e qp 



Euclidean ; 
Lorentzian . 



(92) 



We would like to emphasize that Cq is the sound speed in 
region 2 for k — > 0. 

Again, we take the previous results found in the 
hydrodynamic limit, see equations (|69l l70l ITTj) . where 
we need only replace X-^- with X*, and lo 2 - with uj 2 -. 

Altogether we expect the mean number of particles 
produced for \k\ > c^/e qp to be equivalent to those from 
the purely Lorentzian variations, while phononic modes, 
with \k\ "C c^/eqp, should experience exponential growth 
similar to the result from the hydrodynamic limit. We 
have plotted the mean number of particles produced for 
the whole event, see figure [4(b)] The plots are m agree- 
ment with our predictions. 

One motivation for including microscopic corrections 
into our hydrodynamic calculations was to solve the prob- 
lem of an infinite number of particles being produced in 
sudden variations; see Eq. ([80]) . To estimate the total 
number of particles, we need only to focus on the region 
\k\ > qj/e qp . It can be seen that for large k the mean 
number of particles scales as Nk ~ fc~ 4 . Therefore, for an 
isotropic emergent spacetime (with two or three spatial 
dimensions d), we get 



<k 



Z 7T 



dk k 



d-5 



k d - 



(93) 



(91) 



Therefore, the total number of particles produced is fi- 
nite, as expected. However, it is easy to see that the total 
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energy emitted, 



E = 2 d -\ 



I IT 



dk k 



d-3 



(94) 



is still infinite. Note, that in the ultraviolet regime one 
has to use the non- linear dispersion relation uj k ~ k 2 ; see 
Eq. ([55)1 . We will discuss this further in the conclusion, 
after we investigate a possible application of signature 
change events for laboratory cosmology. 



VI. PARTICLE AMPLIFIER FOR 
COSMOLOGICAL PARTICLE PRODUCTION 

In a Euclidean-signature emergent spacetimc some 
modes (depending on the strength of the attractive inter- 
actions) grow and decay exponentially. As pointed out 
in [2(| this behavior is in analogy to cosmological parti- 
cle production, where super-Hubble horizon modes (i.e., 
modes with frequencies that are smaller than the Hub- 
ble frequency) show similar kinematics. These modes are 
not free to oscillate, as they get dragged along with the 
spacetime fabric. This motivated us to investigate how 
pre-existing condensate perturbations are influenced by 
the existence of a short-duration Euclidean phase. To be 
more specific, we would like to determine whether sig- 
nificant amplification of a pre-existing particle spectrum 
occurs from exposing it to a short finite time of attrac- 
tive interactions. This might be of interest for laboratory 
cosmology in Bose gases. There one of the outstanding 
problems involved with such experiments, as suggested 
in [SI, [33, 111 [Si, SO, S3, is the detectability of the 
quasi-particle spectrum in an emergent FRW-type uni- 
verse. Due to the smallness of mode-population, common 
detection mechanisms would fail to measure the spectra. 
Recently there has been some attempt to solve this prob- 
lem involving a more sophisticated detection mechanism 
[42! ]. Another way to tackle this problem is to amplify 
the quasi-particle spectrum so that common detection 
mechanisms can be applied. 



A. Main concept 

In Fig. ([3]) we illustrate the principle of our particle am- 
plifier idea. The particle amplifier involves a 3-step pro- 
cess. Region 2 will create some quasi-particle spectrum 
fa = /3 IN , which we wish to amplify. The Euclidean re- 
gion 3 is the core of our particle amplifier process with 
fa = /3 PA , and region 4 contains the amplified quasi- 
particle spectrum, fa = /3 OUT . For simplicity we focus 
on a Lorcntzian sudden step process as the source of our 
initial quasi-particle spectrum, however we stress that 
this method is more general and applies to any quasi- 
particle spectrum in region 2 that is of a FRW-type uni- 
verse. 

We would like to determine an expression for our final 
amplified spectrum, in terms of the input spectrum, and 
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FIG. 3: The figure shows a particle production process start- 
ing from the vacuum in region 1. The bold line corresponds 
to particles produced in a single Lorentzian step from region 
1 — » 2, while the dotted line shows particles produced during 
region 2. Region 3 depicts a subsequent amplification process 
of a finite time Euclidean interval. Altogether the OUT signal 
in region 4 can be viewed as an amplified IN signal. 



the spectrum we would expect from the amplifying step. 
We know the connection matrix after the amplifying step 
can be decomposed as (see Eq. dBU)): 



M— = M— ■ M— 

24 12 . ^24* ^12 12 ^24 . n 12 24* 

a k a k + Pk Pk a k Pk + Pk a k 

0I2* 24 . 12* ,o24* i0 12* /3 24 . 12* 24* 

Pk a k + a k Pk Pk Pk + a k a k 



.(95) 



We now look at an expression for the final output spec- 
trum in terms of the input spectrum and particle am- 
plifier contribution. As we already have the Bogoli- 
ubov coefficients, this is straightforward to calculate. As 
det(M*) = 1 it follows that: 



l/3f| 2 



\€ | 2 + [2|/3f | 2 



l]|/?f| 2 



n n I 24 * a 12 a 24 -* 12 * 

2Rc [ajr fa fa ajr 



(96) 



This can be rewritten completely in terms of IN, OUT 
and PA formalism, 



\0 



OUT 1 2 



3PA12 



= |/3 PA | 2 + [2|^ 
+ 2sin( 7 ) |a PA 



1 



3lN|2 



.IN 1 1 0IN1 



(97) 



where 



sin( 7 ) = Re [ph{ar)ph{fa)ph{fa~)ph{aj) , (98) 



12 



represents the relationship between the phases, ph(...), 
of our Bogoliubov coefficients. Therefore the amplifica- 
tion process depends on 7, and it is necessary to find an 
explicit expression for it. 

Although it would be possible to make further state- 
ments while keeping the analysis completely general, we 
now choose an explicit particle amplifier configuration 
and provide an example that demonstrates the parame- 
ters 7 depends on. 

B. Simple example 

For the simplest case, we look at a Euclidean ampli- 
fying step, with X— = i, X M = — i, X— = 1 and write 
= This describes the scenario where the speeds 

of sound before and after the Euclidean step are equiv- 
alent. With these replacements the Bogoliubov coeffi- 
cients for the Euclidean amplifying step can be split up 
as follows: 

a M = ei ^ + T±)-i4(t^ + T±+Tl) |q Mj ^ (g9) 

p m = e+i( 4 ( tia + ^)+i4(^+ri + T2) |/3 M| _ (100) 

Using the notation for the Bogoliubov coefficients de- 
scribing the input spectrum step to be; 

af = \af\ , (101) 

Ht = \/3f\ , (102) 

we then find the output spectrum after amplification be- 
comes: 

\@t\ 2 = i/?r i 2 + (21c i 2 + 1) m 

-2 sin (2w| T L ) (103) 

x v%fl 2 + i)(l/3fl 2 + i) l/fl K\ > 

where we have employed the normalization; det(M 4i ) = 
1. In this particular case we found 7 = —2uj~^T k , which 
therefore only depends on known parameters — the time- 
intervall T 1 and the dispersion relation lo\ in region 2 — 
that are tunable (for a certain range of fe-values) through- 
out a BEC experiment. 

Depending on 7 the number of OUT particles oscillates 
around 

l/3 OU T = (l + 2|/3 IN | 2 |/3 PA | 2 ,(104) 

and thus the difference between a process with and with- 
out initial particles is given by, 

|/3 OUT | 2 a 2|/? IN | 2 |/3 PA | 2 . (105) 

In the last step we assumed that |/3 IN | <§C |/3 PA |, which is 
plausible since we expect a much larger particle produc- 
tion in the Euclidean region. 



Nevertheless, we have shown that for the simplest ex- 
ample possible the particle amplifier does indeed amplify 
the initial spectrum, see Eq. (| 105|) . However, consider- 
ably more effort is required to analyse the particle ampli- 
fier in depth. It will be important to explore how a more 
realistic initial particle spectrum resulting from cosmo- 
logical particle production is amplified and the efficiency 
of our proposed model. The difference between the max- 
imum possible particle production from the amplifying 
step (sin(7) = — 1) and the minimum particles produced 
(sin(7) = 1) might be a possible way to increase the ef- 
ficiency of the particle amplifier. In the next section we 
will review our results, connect with the existing litera- 
ture, and discuss problems and open questions. 

VII. CONCLUSIONS 

In the first few pages of this article we introduced the 
reader to the concept of emergent spacctimcs. Under cer- 
tain conditions, the microscopic theory of an ultra-cold 
weakly interacting Bose gas gives way to a mean-field 
description, the Bose-Einstein condensate. A new state 
of matter has then been formed, described by a com- 
plex classical field. Its behavior is dominated by collec- 
tive variables, rather then individual atoms. A careful 
study of the quantum excitations around the mean-field 
reveals the correspondence between the condensed mat- 
ter physics and quantum field theory in curved space- 
times. The quantum perturbations are massless spin-zero 
particles, that are controlled by an emergent geometri- 
cal/gravitational field. 

The actual excitations 5ijj and 6%/)' are related to 
perturbations in the collective Hermitian variables, the 
phase <f> and density h of the complex mean-field. The 
condensed matter description for linearized perturbations 
supplies commutation relations for the phase and (up to a 
function) its material derivatives. We were able to show 
that these commutators are precisely the ones for the 
field 4> and its canonical momentum II; on the emergent 
curved spacetime. 

The effective gravitational field defined on the Bose 
gas naturally implements a mechanism to experimen- 
tally perform signature-change events. The signature of 
the emergent gravitational field is directly correlated to 
the sign of the s-wave scattering amplitude of the micro- 
scopic collisions. From a gas with repulsive/ attractive 
interaction we expect an emergent Lorentzian/ Rieman- 
nian geometry. The existence of Riemannian geometries 
in a Bose gas is limited to short duration time intervals 
and weak interactions, since back-reaction and instabil- 
ity issues would otherwise destroy the condensate; see 

Inspired by recent condensed matter experiments 1 1 71 . 
[HI , and theoretical work based on them 0, HO, Hi], |43j . 
we investigated the quasi-particle production caused by 
an effective gravitational field going through a sudden 
but finite duration transition to Euclidean signature. We 
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(a) Hydrodynamic limit; L-L-L. (b) Microscopic corrections; L-L-L. 




(c) Hydrodynamic limit; L-E-L. (d) Microscopic corrections; L-E-L. 



FIG. 4: The figure displays the mean number of particles Nk produced from two sudden variations in the strength of the 
atomic interactions. The different curves (color online) correspond to the time period T for which the system is kept in the 
intermediate regime. The left column focusses on purely collective effects, while in the right column microscopic corrections 
are taken into account. The top two diagrams show variations between spacetimes emerging from a Bose gas with repulsive 
interactions (a > 0). In the bottom row of the diagrams the intermediate regime exhibits an Euclidean geometry, where the 
underlying Bose gas shows attractive atomic interactions (a < 0). 



compared the signature changing case with sudden vari- 
ation between different Lorentzian regimes, and showed 
that perturbations in the Euclidean regime experience 
exponential growth, while for an intermediate Lorentzian 
regime the modes oscillate between single steps. These 
are given by the maxima and minima configuration of 
the two-step process with T = 0. Our results are in 
two ways a generalization of the existing literature, (e.g., 
the standard general relativity calculations of Dray et al. 

mm): 

1) For the right choice of mode functions (a complex 
conjugate pair), the calculation is formally equivalent 
to purely Lorentzian sudden variation between arbitrary 
levels. This formalism can be applied to any n-step pro- 
cess. 

2) In addition, we went beyond the hydrodynamic de- 



scription, which is only sufficient to account for the in- 
frared behavior of the system, and included ultraviolet 
modes in our analysis. This leads to modifications in 
both the emergent gravitational field and the dispersion 
relation. For the gravitational field we obtained (in the 
eikonal limit) a momentum-dependent rainbow geome- 
try; and a non-relativistic dispersion relation. The mod- 
ification in the dispersion relation can be classified as be- 
ing within the boost sub-group, and is of the form as sug- 
gested by some effective field theories [H H, [H [H, H| . 

In figure (0]) we illustrated our four different results, 
namely sudden variations without (first row) or with (sec- 
ond row) signature change events, within the hydrody- 
namic limit (left column), and beyond it (right column). 
As a result, the total number of particles produced is 
finite. One way to explain the behavior of the trans- 
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Planckian modes is to understand the physical signifi- 
cance of the healing length of the Bose-Einstein conden- 
sate. The healing scale defines a length scale over which 
quantum fluctuations in the condensate tend to smooth 
out. In some sense, it defines the smallest possible size 
for a condensate, V > see (35J. Perturbations with 
wavelengths that arc much larger than the healing length 
experience a nice smooth mean-field, here our emergent 
gravitational field. High energy modes, with wavelength 
at the magnitudes of the healing length start to see fin- 
gerprints of the microscopic structure. Note that, per- 
turbations of sufficiently high energy are not driven by 
collective variables. This can also be seen in figure [4(d)] 
While infrared modes behave like their hydrodynamic/ 
relativistic neighbours to the left 4(c) the ultraviolet 
modes do not notice the change in the signature of the 
gravitational field. They arc behaving just like their up- 
per neighbours |4(b)| in a purely Lorentzian geometry. 

Still, not all of our concerns are yet resolved. For ex- 
ample, while the total number of particles produced can 
be made finite, the total energy remains formally infi- 
nite even beyond the hydrodynamic limit. This problem 
might be related to our choice of ultra-high-energy de- 
scription for the fundamental Bosons. Perturbations of 
sufficiently high energy can excite single atoms out of 
the condensate, and these are (in our description) non- 
relativistic particles. Therefore, we suggest that to im- 
prove the ultra-high-energy behaviour it might be useful 
to start with a fully relativistic description for the fun- 
damental Bosons. (For a description of relativistic Bose- 
Einstein condensates see, e.g., |44|.) 

Finally, we would like to comment on the possible con- 
nection between our theoretical results and the data from 
the Bose-nova experiment carried out by Donley et al. 
[l7l . [l8| . There are fundamental differences between our 
approach, and the Bose-nova experiment. First, we have 
chosen a hard-walled box as an external potential, while 
in the experiment a harmonic trap has been used. In 
@, US EL E3], Hu and Calzetta have shown that the 
trapping energy delays the condensate collapse for a cer- 
tain amount of time. Their calculations are in good ac- 
cordance with the experimental values. In our case the 
trapping energy is zero, and hence we do not expect such 
a delay. 

Another fundamental difference is that we kept our 
gravitational field non-degenerate (apart from the actual 
instant of signature change) . In the experiment the initial 
scattering length was taken to be zero and held there for a 
finite time, and hence Donley et al. started with a partly 
degenerate gravitational field. (Partly degenerate gravi- 
tational field, because for a harmonic trap the quantum 
pressure term exhibits a fc-independent contribution re- 
sulting in a position-dependent correction to the effective 
atomic interaction; see first term inside the curly brack- 
ets in Eq. I|82p.) In the experiment they have chosen two 
different sequences for the sudden variations, with two 
different outcomes. One, where the scattering length has 
been driven from zero to attractive and finally to large 



repulsive atomic interactions. The other set-up was from 
zero to attractive, back to zero and then to large re- 
pulsive atomic interactions. In both cases they detected 
bursts of atoms leaving the condensate. But the bursts 
without the intermediate regime of zero interactions were 
much stronger. These are "jets" of atoms, while in the 
other case the bursts were less strong. Hu and Calzetta 
approached this problem without taking back-reaction 
effects — of the quantum perturbations onto the con- 
densate — into account. Naturally, their description 
becomes less accurate with the length of the Euclidean 
time-interval. We would like to emphasize the difference 
in the behavior for the acoustic metric, the "gravitational 
field" , since in the first scenario it does not exhibit any 
finite interval of degeneracy, while in the latter it does. 

Furthermore, we would like to propose a rather dif- 
ferent use for our L-E-L-process, as a particle ampli- 
fier for cosmological particle production in the labora- 
tory. In our calculations so far we always started from 
the vacuum with zero collective excitations in the sys- 
tem. These calculations could easily be extended, such 
as to start with a non-vacuum state, e.g., after cosmo- 
logical particle production in an emergent Friedmann- 
Robertson-Walker type geometry [H, IS3> HI S E El • 
Recently this behavior has been studied numerically in 
a realistic Bose-Einstein condensate, where the present 
authors used a classical-field-method approach, see [28| . 
To obtain an effective expanding universe, the scattering 
length has to decrease as a specific function of time. Af- 
ter a finite expansion time, the quasi-particle spectrum 
obtained — a very small number of particles — somehow 
has to be detected. Due to its smallness this remains a 
significant experimental challenge. Recently, there has 
been some theoretical effort regarding this problem, see 
[4^ . That author suggests a rather complicated detec- 
tion mechanism, where only perturbations of one wave- 
length at a time can be detected. In section (|VI[) we 
suggested an alternative concept, that is to amplify the 
quasi-particle spectrum, e.g., in a three-step process in- 
volving a brief finite-duration Euclidean region to am- 
plify the signal, so that common detection mechanisms 
might be able to detect the amplified quasi-particle spec- 
trum. For a simple example, we calculated the particle 
spectrum after the amplification process in terms of the 
initial spectrum. However, for any practical application 
in a real BEC experiment there are considerable techni- 
cal issues that need to be explored. (For example, the 
possibility of using a L-L-L amplifier.) 

Before finishing we would like to briefly discuss the pos- 
sibility of real general relativistic signature change events. 
In 1992, James Hartle and Steven Hawking proposed a 
signature change at extremely early times, in the very 
early stages of the big bang, when quantum gravity ef- 
fects are expected to be dominant (see reference |45j). 
They suggested that the existence of physical time, and 
hence the existence of our universe, is associated with 
a signature change event from Euclidean to Lorentzian 
geometry (the "no boundary" proposal). 
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In some sense it is possible to consider the reverse Bose- 
Nova experiment (attractive to repulsive) as the creation 
process for our emergent spacetime. At early times the 
Boson interactions are strongly attractive, such that the 
atoms do not show any collective/ mean-field behavior. If 
now the interactions experience a smooth, or non-smooth 
change for the atomic interactions, from attractive to re- 
pulsive, the individual atoms have to give way to collec- 
tive/ mean-field variables. Therefore the existence of our 
emergent gravitational field might in some sense be asso- 
ciated with a pre-dating signature change event caused 
by a change in the underlying microscopic variables. 
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APPENDIX A: JUNCTION CONDITIONS 

For the benefit of readers who may not wish to deal 
with differential forms and conjugate momenta, we now 
present a "low-brow" calculation of the a and (3 coeffi- 
cients — for the specific physics problem we are inter- 
ested in — in two different coordinate systems. Both 
co-ordinate systems lead to the same physics result, and 
both give explicit expressions for the conserved inner 
product, the modes, how to normalize them, and of 
course the junction conditions. 



subsection where we explain how we could have equally 
well done things in terms of "canonical time" . One can 
read the conserved inner product off by starting directly 
from the PDE itself — noting that this PDE is defined 
for all time, and then considering the quantity 

{KAb) = j (^-jj -fajj ddx > ( A4 ) 

where <f> a and </>& are solutions of the master PDE. Then 
by Gauss' law (or the fundamental theorem of calculus) 



Oa-, 9b ) I final 

dt 



l>a: 9b ) initial 



(A5) 



<t> a jj d t <t>b - <t>b jj d t (j)^j d d+1 x. 



Now apply Leibnitz' rule, the PDE, and an integration 
by parts to obtain 



t>a, <Mfinal ~ (<j>a, 4>b)\ 



itial 



0. 



(A6) 



So this is the correct conserved inner product for the 
Klein-Gordon-like PDE we are starting with. Note that 
this conservation law, the way we have set it up holds 
for any and all initial and final times, regardless of which 
side of the junction they are located. 



2. The two non-overlapping simple regions 

In region 1 (t < t—) the master PDE reduces to 



^ V 2 = 0, 
m 



(A7) 



Basics 



The master PDE that everything starts from is 
\ U I m 



(Al) 



where n , m and h are constant, while 

U(t) = U 1 -Q(t^-t) + U 1 <d(t-t^). (A2) 

The most basic forms of the junction conditions are then 
that the field operators </>, and its conjugate momentum 
(h/U) d t (j>, are continuous at t—. If these conditions are 
not satisfied, then there is no way that the master PDE 
can be satisfied at t— — there would be delta-function 
contributions on the LHS that would not cancel against 
anything on the RHS. These two junction conditions can 
be written as 



[<t>\ = 0; 



— Ot(, 

U 



= 0. 



(A3) 



Note that we are here working in terms of physical labo- 
ratory time, and will continue to do so until we get to the 



and can be rearranged to 



m 



Introducing 



this becomes 



i\2 



n U± 



dfy + (c 1 ) 2 V 2 = 0, 



(A8) 



(A9) 



(A10) 



with the understanding that this PDE holds only in re- 
gion 1 (i.e., t < t—). Similarly, in region 2 (i.e., t > t—) 
the master PDE reduces to 



d^-Kc 2 ) 2 V 2 ^ = 0. 



(All) 



The junction conditions are unchanged, though for con- 
venience we can write them in terms of c 2 = no U/m 
as 



[0] = 0; 



0. 



(A12) 
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In each individual region the field satisfies the usual 
flat spacetime Klein-Gordon equation, so solutions are of 
the form 

(p oc exp(i[uJkt — kx]) , (A13) 

where u>k and k satisfy the dispersion relation 

^- = cfc, (A14) 

and we have very carefully not yet specified any normal- 
ization for these modes. (Nor is c the same in the two 
regions, it is either c- or c- as appropriate.) 

One thing we can say without further calculation is 
this: In view of the junction condition [<f>] = 0, that is 
(j)(t— ~,x) = <j){t— + ,x), the spatial position dependence 
of the solutions to the PDE on the two sides of the junc- 
tion must be the same — this implies that in terms of 
the plane waves above we must enforce k to be the same 
in regions 1 and 2, and therefore 



Normalization 



(A15) 



Now let us compare the conserved inner product for the 
PDE we are physically interested in with the inner prod- 
uct for the naive Klein-Gordon equation. In particular, 
if initial is before the transition and final is after the tran- 
sition, then in terms of the naive ordinary Klein-Gordon 
inner product the conserved inner product for our PDE 
is 



K , 4>b ) I initial = JJl {(j>a,<i>b) 



(4>a, <M|final 



where (<f> a , </>&) na i ve = / {(pa dt<t>b - (f>b <9 t Q ) d d x . 
This means that properly normalized modes are 



U 1 



cxp(i[uj k t - kx]), 



(A16) 
(A17) 

(A18) 



which (ignoring a trivial overall constant factor, that does 
not change from region 1 to region 2) we might as well 
write as 

exp(i[uj k t - kx]). (A19) 



V2w/T 

In particular in regions 1 and 2 we want to deal with 

c 



and 



vr k (t,x) 
respectively. 



exp(i[w7i-kx]), (A20) 



exp(i[w|i-kx]), (A21) 



4. Applying the junction conditions 

In regions 1 and 2 we write the solutions of the PDE 
as the real parts of 



A 



+B- 



C7. 



exp(i[w^ - kx]) 



exp(i[-u^-kx]), (A22) 



and 



C- 



+D 



exp(i[w|t — kx]) 



exp(i[-w|t- kx]). (A23) 



Note that the x dependence is the same in all four of 
these terms and so quietly factors out — this is why we 
asserted that k had to be the same in both regions. 

Applying the first junction condition [0] = 0, at the 
transition time t— , and using the dispersion relation im- 
plies 

Vc 1 [A exp(iw^i il ) + B expi-iujlt 11 )] 

' [C exp{iur k t m ) + D exp(-iw|^)] . (A24) 



The second junction condition, [9 t 



0, applied 



at the transition time t— , implies (after factoring out all 
the c's and uj's, and using the dispersion relation) that 

—= [A cxp(ia^£— ) — B exp(— iuj\t— )] 

[C expiiult^) - D exp(-iuj g k t^)] . (A25) 



These two junction conditions are enough to com- 
pletely specify the transmission matrix of a's and /?'s. 
Indeed, setting A — > 1, B — > 0, C — > a, and D — > /3 and 
solving for a and /? we find 



cxp(-i[c 2 -c 1 ]H i2 ); (A26) 




13 =o \ \ - -\ - exp^ + c^); 



(A27) 



which is exactly the result reported in the body of the 
paper. 

5. Using "canonical time" 



Let us now pick another time coordinate. Let us call 
it "canonical time" and define it by 

T=l C Mdt= f^dt^ (A28) 
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= t-Mtt dt = S ^ dt, 



(c.) s 



(A29) 



where we have introduced a convenient constant reference 
point U* and used this to define a convenient constant 
reference speed 



2 _ n [/* 



(A30) 



The the master PDE, which was in the original labo- 
ratory time coordinate 



U J m 



now becomes 



(A31) 



(A32) 



Recognizing that noU 2 /(m U) = (c») 4 /c 2 and defining 

(c*) 2 



C e ff(T) 



c(T)' 



(A33) 



this can be re-written simply in the form of a "paramet- 
rically excited oscillator" 



d 2 (f> + c 2 cS (T) V 2 = O. 



(A34) 



Assuming an exponential space dependence, separa- 
tion of variables yields 



so that 



4>(t,x) = <f>{t) exp(— ikx), 



dU = c 2 cS (T) k 2 I 



(A35) 



(A36) 



The junction conditions, which were originally in terms 
of laboratory time (refer to equation |A12| ) now become, 
in terms of canonical time, the very simple: 

[0] = 0; [d T <t>] = 0. (A37) 

The conserved inner product (|A4[) can now be re-written 

J- n 

(<f> a: 4> b ) = — j (0 a d T <f> b - <j> b d T K) d d x, (A38) 



that is 



t>a, 0b) = jj- {<t>a, <pb)t 



(A39) 



So in "canonical time" coordinates we have a particularly 
simple PDE, elementary junction conditions, and a triv- 
ial inner product — the only "tricky" thing is that we 
have to use the "effective" sound speed c c g = (c») 2 /c. 

Now in region 1 the original master PDE can be cast 
into the form 



-a 



(c*) 4 



V 2 = O, 



(A40) 



while in region 2 



(c*) 4 

( C 2) 2 



V z (h = 0. 



(A41) 



In view of the form of the conserved inner product, the 
normalized modes are 



1 



exp(i[wiT - kx]); 



exp(i[^T-kx]); 



(A42) 



(A43) 



Note the "odd looking" form of the dispersion relation 
- but this is just because ui^ and w| are not physical 
frequencies — they are "T-time frequencies" . 

In regions 1 and 2 we write the solutions of the PDE 
as the real parts of 

exp(£[c^T - kx]) 



-B 



= C- 



+D- 



exp(i[-wtT - kx]); (A44) 



cxp(i[uj 2 k T — kx]) 



exp(i[-w^r-kx]). (A45) 



The first junction condition [0] = 0, after applying the 
dispersion relation is 

[A exp^T 12 ) + B expH4T^)] (A46) 
[C exp(iw£r-) + D exp(-^|T i1 )] . 

The second junction condition [dr^] = 0; after simpli- 
fying using the dispersion relation, now leads to 

-= [A exp(«jiT a ) - B expHc^T 12 )] (A47) 
yc 1 

= -$=[C exp(ic^T-) - D exp(-z^T^)] . 
Vc 2 



Up to irrelevant phases, these are the same equations as 
were derived in laboratory time. 

Now setting A — ► 1, B — ► 0, C — ► a, and D — ► /? and 
the solving for a and /3 we find 



1 



+ J- exp^-u^), (A48) 



exp(i[w£ +w|]T^). (A49) 
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The phases on a and (3 are now slightly different, but 
this is not new physics — it has to do with the phases we 
picked for our "normalized modes" — these phases are 
now slightly different from the laboratory time calcula- 
tion. Note the magnitudes are completely unambiguous: 




(A50) 




(A51) 



In short, this appendix has served to verify that the key 
technical parts of the calculation can mathematically be 
carried out in a number of different ways that ultimately 
lead to the same physical result. 
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